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Abstract 

A time dependent variational approach is considered to derive the equations of movement for the 
A(/)^ model. The temporal evolution of the model is performed numerically in the frame of the Gaussian 
approximation in a lattice of 1+1 dimensions given non homogeneous initial conditions (like bubbles) for 
the classical and quantum parts of the field which expands. A schematic model for the initial conditions 
is presented considering the model at finite fermionic density. The non zero fermionic density may lead 
either to the restoration of the symmetry or to an even more asymmetric phase. Both kinds of situations 
are considered as initial conditions and the eventual differences in early time dynamics are discussed. 
In the early time evolution there is strong energy exchange between the classical and quantum parts of 
the field as the initial configuration expands. The contribution of the quantum fluctuations is discussed 
especially in the strong coupling constant limit. The continuum limit is analyzed. 
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1 INTRODUCTION 



There are many motivations for the study of time dependent non perturbative methods in Quantum Field 
Theory (QFT). Some examples are immediately found in systems with strong coupling constant, sponta- 
neously symmetry breaking and which undergo phase transitions. Some of the most interesting cases are 
present in the relativistic heavy ion collisions currently being prepared and performed in BNL/RHIC which 
probe hadronic matter at very high densities and temperatures. In these cases a region of very high energy 
density starts expanding (and "cooling") just after the nuclei collision. These systems are usually described 
by hydro-dynamical models which are known to be quite reliable [|^]. However, the understanding of these 
models in terms of microscopical descriptions with the underlying QFT degrees of freedom are expected 
to be derived and are currently being investigated by several groups. Many effects are expected to occur 
in those systems such as, for example, dynamical phase transitions whose existence should be present in a 
reliable description. 

Due to the extreme complexity of realistic theories, as QCD, one usually is lead to study effective models 
which respect the major properties of the fundamental theory and reproduce the main issues of it in some 
range of a relevant variable (as energy). In the present work, however, a still more simplified version of 
the reality is considered in order to check qualitative effects. The model is often used as a test model, 
although its scalar field may be considered as the relevant degree of freedom for inflationary models in 
Cosmology ||2|. In condensed matter and statistical mechanics it is also of interest [^. Furthermore it can 
be identified to the mesonic sector of the linear sigma 0{N) model in the large N limit or without pions. 

In the last decade a quite large amount of work have been done in order to shed light in some of 
the subjects mentioned above as well as aspects related to the dissipation, thermalization, decoherence, 
formation of disoriented chiral condensates (DCC), phase transitions among others with non perturbative 
time dependent formalisms. Some examples are found in [§, |, |, 0, |8|, |^], ||ll|, |l2|, |l5|, Some 

works have already been done concerning the dynamics of non homogeneous configuration in bosonic fields. 
In particular, the Gaussian equations of movement were considered in ||8|, ^, 18] to study properties of the 



large time dynamics in certain cases. Non homogeneous field theory was also studied in |18| to obtain 
information relevant for particle production and thermalization of non-equilibrium systems. 

In the present work we extend, complement and give a sounder picture for the work done in [^] . The main 
aspects we study in the present work are the following. Firstly a schematic model for the formation of locally 
non equilibrium initial conditions (which expands with time) is developed. As already proposed in this last 
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paper, we consider the possibility that, due to some particular condition, the condensate amplitude (as well 
as the physical mass of the scalars) is either suppressed or enhanced in a localized region of space. For this 
we consider the Xcj)^ model coupled to fermions at finite density in a one dimensional space. Depending on 
the kind of the coupling the model may experience either symmetry restoration at high fermionic density or 



further symmetry breakdowns |16]. We look for dynamical consequences of the corresponding enhancement 
or suppression of the condensate at the tree level and in the frame of the Gaussian approximation. Although 
there are enormous differences between this simple (idealized) model (the Xcj)'^ model) and the realistic 
"fireballs" from RHIC, we believe that the study of the influence of the present field theoretical degrees of 



freedom is of interest. A partially similar idea to this one was discussed in reference [17| where the pressure 
due to a gas of pions, including a X(j)'^ term in a static description, was considered to drive an expansion of 
a "fireball". 

We have therefore the following picture. Firstly we fix two parameters of the model, which will allow 
us to make meaningful comparisons between the (time dependent) tree and Gaussian levels. Secondly, it is 
assumed that the scalar field is locally placed in a thermal bath and/or experiences an interaction with a 
finite fermion density in that small region. These interactions -which change the ground state of the model- 
are switched off at t=0 yielding non equilibrium initial conditions for the scalar field which are evolved in 
time. The temporal evolution is performed within the tree level and Gaussian approach equations, producing 
the expansion of regions (bubbles) endowed with high energy densities. 

The work is organized as follows. In section 2 a time dependent variational method for pure states systems 
is described in the Schrodinger picture and the equations of movement are derived with a Gaussian trial 
wave functional. Some considerations for the static and thermal case are done. The small amplitude motion 
case is investigated for homogeneous configurations in order to provide some useful results for analyzing the 
relevance of the one dimensional lattice simulations for a more realistic situation. In section 3 the numerical 
method using a pure states generalized density matrix - Time Dependent Hartree Bobogliubov (TDHB), 
as developed in is extended for the asymmetric case (0 7^ 0) in a discretized space. It allows for the 
investigation of the temporal evolution of a localized region of the space (lattice) where a high energy density 
occurs. Next, in section 4, by coupling the scalar field to a finite fermionic density system we are able to 
construct a model for the initial conditions at finite density. The possibility of symmetry restoration and 
further asymmetric phase(s) is discussed and the temporal evolution of corresponding configurations will 
be analysed. Still in the section 4 we provide an alternative way of fixing the parameters of the model in 
order to perform plausible comparisons between classical and quantum field systems. The numerical results 
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for the early time evolution are presented in section 5 for different initial conditions and values of the free 
parameters of the model. All the numerical examples shown in this article conserve total energy. Special 
attention is given to the strong coupling limit and one example of the temporal evolution of a deviation 
from a Kink solution is exhibited. The continuum limit is discussed. The results are summarized in the last 
section. 



2 GAUSSIAN APPROXIMATION 

The time dependent variational approximation at the Gaussian level has been studied for several years 
1^, 20 1 . It provides a systematic method to study the temporal evolution of a quantum field 

theoretical model given an initial condition by means of the equations of movement. In spite of recent 

achievements for considering approximations beyond the Gaussian (for example in [11, ^) we want to 

address some unexplored features in the time dependent Gaussian approach. 

Let us consider the time-dependent Dirac's variational principle with a trial Gaussian wave functional 
m. Firstly the average value of an action / is calculated with a given trial wave functional |^ >: 



I = J dt <'i'\{idt- H)]"^ >, (1) 

where H is the Hamiltonian. 

In the Gaussian approximation at zero temperature the wavefunctional ^ is parametrized by: 

^ [0(x)] = iVexp|-^ J dxdy<5</)(x) (g-1(x, y) + iS(x, y)) 6<P{y) + i J dxvf (x),5</.(x)| . (2) 

Where 5(j){x,t) = <^(x) — (j){x,t); the normalization is A^, the variational parameters are the condensate 
(j){x,t) =< ^'101^' > and its conjugated variable 7f(x, t) =< ^'|7r|^' >; quantum fluctuations represented by 
the width of the Gaussian G(x, y,t) =< ^\(f)(x.)(f)(y)\^ > and its conjugated variable $](x, y, t). 

The Lagrangian for a scalar field cp with bare mass rrig and quartic coupling constant A is given by: 

£(x) = i |5^<A(x)9^<A(x) - mlH^f - ^<P{^)'] ■ (3) 

Prom this expression the corresponding Hamiltonian H is obtained. The action of the field operator (p and 
its conjugated momentum vr in functional Schroedinger representation over a wave functional ^ [4>i^)] =< 
0(x)|^ [(j)] > is respectively: 

<^|^[<^(x)]>=0(x)|M/[0(x)]>, 

7r|^ [0(x)] >= -i(5/(50(x)|^ [(/.(x)] > . 
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The average value of the Hamiltonian, which will be explored in the numerical simulations, is given in terms 
of the variational parameters with the functional integrations: 

iG-i(x,x) + 4SGS(x,x) + 7f2(x) - AG(x,x) + m2G(x,x) + ^G2(x,x)+ 



n = ^ 



+ (V0(x))^ + A^4(^) + t0^(x)G(x, x) . 
In the Schrodinger picture the wave functional evolves like the Schroedinger equation: 



(5) 



i^*[0(x)]=F^[0(x)]. (6) 

This equation is equivalent to the temporal evolution of the variational parameters given initial conditions. 

The variations of the / with relation to the variational parameters (of the Gaussian wave functional) 
produce the equations of movement. 

61 



5E(x,y,t) 



dtG{x, y,t) = 2 (G(x, z, t)S(z, y, t) + S(x, z, t)G(z, y, t)) 
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atS(x,y,t) = (^2S(x,z,t)S(z,y,f) - -G-^(x,y,i)j + 



a,0(x,t) = -vf(x,t) 
atvf (x, t) = r(x, y, t)4>{y, t) + t) 



(x,t) ^ '^^ ^ '-^-^^ 6 

where: r(x, y,t) = —A + (rriQ + ^G(x,x, t)j 6{x — y). In the x-coordinate space the Gaussian width is 



written as: 

G(x,y)=<x|- ^ |y>. (8) 
V-A + 

As it is well known the static limit of this approximation produces the "Cactus" diagrams resummation 
p^]. In section 3, an alternative way of writing these equations will be discussed and evolved in time in 
a lattice for a class of non homogeneous initial conditions. It corresponds to the time dependent Hartree 
Bogoliubov approximation 

By performing a Fourier transformation it is possible to eliminate the variables vf and E. The equations 
in the asymmetric phase become: 



Gkk' - - +2(k^ + ml + ^G{x, x) + ^0 ) G^t 

k + (fc^ + ml + ^<^2 + ^_G{x, x)\ 4)k = 



(9) 



with Gjfj^/ =< k|G(x, x)|k + q >. These equations were generalized for the out of equilibrium (non zero 
temperature) using different methods in [Q, 10, 19, 21 1 and they were studied extensively mainly by means 



of numerical calculations. The choice of initial conditions is entirely subordinate to the approximation, in 
the sense that were it not Gaussian one might have to consider three conditions instead of two [^]. The 
analysis of these equations in show that initial conditions (for homogeneous G and (f) ) are crucial for the 
time interval in which the system evolves towards the minimum and beyond as well as for the speed of 
the field evolution ||^, ^. Some other effects have been addressed as, for example, dissipation via particle 
production, the Landau damping, collisional relaxation at zero temperature as well as the relevance of the 
initial conditions for the dynamics was also investigated in some cases in the early and late time evolution 
I, |, |, H, |2|]. In particular, it was found that the Hartree approximation is well suited for the 



study of early-time dynamics. 

2.1 VACUUM AND RENORMALIZATION 

The state of minimum energy is found from the equations of movement in the static case G = S = </) = 0. 

The equations of motion reduce to the GAP equations which minimize the effective potential. They can 
then be written as: 

2 2 , ^ 

■ = -I 

2 

These equations provide the two phase Xcf)^ model: a symmetric phase (where there is only a zero condensate 



-r + mi + -G{n'))=Q, 

/ (10) 



= 0) and the asymmetric phase where the condensate in non zero in the vacuum. From these equations, 
in the asymmetric phase, we find that: 



,2 



2 ^ (11) 



3/^ 

In spite of written in terms of the bare coupling constant this value can be compared to the tree level value 
In the former case (expression (pT[)) the bare coupling constant is fixed by the value of the 



A 

renormalized coupling, as discussed below. 

For a thermal environment, in d dimensions, we can perform a similar calculation to the above one with 
a density matrix with mixed states. It yields thermal fiuctuations corrections to the two point Green's 
function G(x, x) which is substituted by: 

£t, 2> f ^2 f * ' + ""'1^ 2 ) f dk 1/1 ,,, A (12) 

In the above expression /(k) is the Bose-Einstein occupation number. It is well known that these thermal 
fluctuations for the asymmetric phase may restore the symmetry: the vacuum solution of the condensate 



4>o tends to become zero as temperature increases. At very high temperatures (T >> /x) the integral (12) 
in one spatial dimension can be expanded and the second expression of (p!o|) can be written in one spatial 
dimension as: 



= " I + ^(Zn(/./(47rT) + 7)] + 0{T% (13) 

where 7 = 0.57721 and /u^ is the mass at zero temperature. Therefore the physical mass is reduced at high 
temperatures, fij^ = /x^(T) < fj?. This effect is also present in the condensate (by means of expression (p^)) 
making possible the eventual restoration of the symmetry breaking (at least at this level of approximation) 

The renormalization of the time dependent approach (in continuum spaces) can be done in the same way 
at zero and finite temperatures as well as out-of-thermodynamical equilibrium [p!o| , [191 , pO| . No additional 
ultraviolet divergence is found due to temporal evolution of the system or its departure from zero temperature 
or equilibrium. Renormalization can therefore be performed by absorbing the divergences in the physical 
mass (1 dimensional system) in the vacuum. By re- writing the GAP equations in terms of the renormalized 
quantities {m'j^, Xr) we can find the relationship between the bare and physical quantities. They can be 
written in one dimension as: 

m^(ld) = mo + —In j 
A 



XR{ld) = A 



1 - \ (14) 



1+ ' 



We note that the coupling constant acquires only a finite correction. 

In a discretized space there is a natural regulator (the lattice spacing for which: A = 1/Ax). The 



integrals (12) become a summation and converge. As the limit to the continuum is taken the integrals tend 



to diverge and a redefinition of the bare quantities is needed. This will be discussed in section 5. 
2.2 SMALL DEVIATIONS AROUND THE MINIMUM 

Before exhibiting numerical solutions of the equation of movement in a lattice let us perform an exercise 
which is useful for the understanding of results. For a certain class of initial conditions it is possible to find 
analytical solutions for the equations of movement. Firstly because we can see the relevance of the temporal 
evolution of a given initial condition irrespectively to the regularization method for the local divergences. 
We can also check that the study of the one dimensional model can provide information about the three 
dimensional case. 
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We therefore consider the symmetric phase of the model {(j) = 0) assmning a particular kind of ho- 
mogeneous initial conditions (exactly the same calculation can be done for non-homogeneous asymmetric 
potential cases, this would not change the conclusions). Let us assume that, for some reason external to the 
model, the system experiences a small deviation from the state of minimum. As discussed before, this can 
happen because of an external thermal bath and it can be parametrized by a change in the physical mass, 
e.g. m{t = 0) = .9/i. This drives quantum fluctuations away from the ground state and generates dynamical 
evolution. In this case it is possible to extract exact analytical solutions for the equation of movement. 
For this kind of initial conditions one can linearize the Gaussian equations of motion using the following 
prescription for the solution: 

where Go(//^) is the value of the fluctuations which obey the GAP equation and 5G{t) the value of the small 
deviation which evolves in time. 

We have presented analogue solutions for the free case in Q choosing a plane wave prescription for the 
deviation 6G{t) since an infinite system is considered. We have found the following solution for the deviation 
6G{t) in d-dimensions: 

where: 

^ J (27r)rf (2Vk2 + ^2)3 3c J ^^^y (2^k2+^2)3 " ^ > 

The integral Is is time independent and contains a log divergence in 3 dimensional space. It can be absorbed 
in the renormalization of the mass and coupling constant. The integral I^c also has a divergence at t = 
in three spatial dimensions which also can be regularized. This problem was also addressed in 1^, In 
particular, the integrals can be written in terms of generalized hypergeometric functions. 

In an one dimensional space, on the other hand, there is no (ultraviolet-UV) divergences for this case 
of temporal evolution. Given a finite initial condition (which is the deviation from the vacuum value, 
proportional to m'^ — ^2) the temporal evolution is finite. This is very significant for the numerical solutions 
in a lattice: the temporal evolution of the system is independent of the regulator. This happens due to the 
fact that the static sector was separated from the dynamic evolution and the divergences only contributed 
for the redefinition of the former. In other words, the regularization and renormalization do not interfere in 
the temporal evolution. This is checked in section (5) with numerical calculations in a 1-dim lattice. 
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Other interesting feature can be pointed out from these solutions by rewriting them. In an 1 dimensional 
space the time dependent integral can be written as: 

_ (fi^-m^) roo j cos(2yt) (.y. 

Whereas in three spatial dimensions the same time dependent integral can be written as: 

r(3) _ (/^^-m^) / 2 roo j^, cos(2yt) nNo(2fit) \ /^^x 
he - 8n [t^ Jf. ^y y2^y2_^2 + 2 ) ' 

The main difference between these two expressions in the fact that in three dimensions there are two 
additional divergences, one of them occurs only at t = 0. This also suggests that the dynamics in one spatial 
dimension is not completely different from the dynamics in three dimensions for the cases under study. 



3 NUMERICAL METHOD: Hartree Bogoliubov 

This section is an extension of the work developed in for the time dependent Hartree Bogoliubov approx- 
imation in the presence of a (time and space dependent) condensate. The generalized Hartree-Bogoliubov 
density matrix can be defined, in a discretized space, as pl|]: 



Ri,j 



(19) 



pi J ^ij 
V ~Pi,j 

where the average density matrices are written in terms of the average of creation and annihilation operators: 
Pi J = ^ < CLiOj + a|aj > which is hermitian and = — < aiUj > is symmetric. 

For the calculation of this matrix, the creation and annihilation operators are written, in a d spatial 
dimensions lattice, in terms of the field and its conjugate variable: 

= -^{mio^)'-^ +i<jm'^} (20) 
at(j) = l={^(j)(a)'^ (21) 

Where a and P are (dimensional) normalization factors. A suitable choice for them is the lattice spacing. 

The averaged values of the field variables ((/>, vr) previously shown can be written in terms of the matrix 
elements of the above generalized matrix R. Using the notation of expressions (|7[^), we obtain at zero 
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temperature: 

< (^{i)(P{j) >= G{i,j) + 1 = JX^^MpiiJ) + ^ihj)) 

< 7r(i)7r(j) >= F{i,j) + 1 = (^^^y+i ^^^ipi^' j) " i^ihj)), ^^2) 

< (A(i)^(i) >= 2G{i,k)^{k,j) +^iTti = -2Im{p{i,j) - 

< 7r{i)^{j) >= 2J:{i,k)G{k,j) +7ti^i = 2Im{p{i,j) + K{i,j)), 

where F{i,j) = G(z, j)^^/4+4S(i, k)G{k,l)Ti{l,i). Therefore we can already note that there are two different 
ways of expressing the state of the system: by means of the generalized density matrix elements OR with the 
variational parameters discussed before. They are completely equivalent descriptions for the prescriptions 
used in this work. 

Next, we discuss the dynamics of the generalized density matrix. The temporal evolution is governed by 
the Hartree- Bogoliubov energy, which can be parametrized in the following form: 

/ 



1 SE 

— Hi 



2' 6Rji 



y -Dij -Wij 



(23) 



It is worth to notice that this approximation, and consequently the dynamics under study, is invariant under 
an unitary transformation. Therefore it is possible to consider the generalized density matrix and Hartree 
Bogoliubov {Hij) energy in another form, namely: 

where the superscript means the transposed matrix and r is given by: 

(25) 



1 -1 
1 1 



With these transformed matrices, one can check that a Liouville-von Neumann type equation is necessarily 
satisfied only in the symmetric phase (i^ = 0): 



Hib, Rk, 



(26) 



This equation can be written in terms of the averaged quantities ( ^2[ ) . By doing this exercise we have shown 
that they are equivalent to equations (0) in the zero temperature limit without the condensate. When 
the classical field is taken into account, the equation ( p6|) is still considered and i^(x) acts as an external 
(dynamical) source to (26). The time dependence of (p{t) is determined by the two last equations of the 
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set (0). These last equations, on their turn, depend on the fluctuations which are evolved in (|2^). We 
are currently extending this method to finite temperature non equilibrium systems [|^]. The numerical 



evolution of (2£) is obtained by diagonalizing Rij and evolving the eigenvectors from which the variational 
parameters (eg. Gij) can be calculated. 



4 PARAMETERS OF THE MODEL AND INITIAL CONDITIONS 

In this section we complete the schematic model for justifying the initial conditions used in the evolution of 
equations of movement. The finite temperature effects were already discussed in the frame of the Gaussian 
approach and now we couple the A(/>^ model to fermions at finite density. There are many possibilities for 
this coupling and we discuss only a few cases which may lead either to the restoration of the spontaneously 
broken symmetry or to a further asymmetric model at high densities. We also provide some remarks which 
are helpful for the choice of the sets of parameters: physical mass fi^ and coupling constant A in the lattice. 

4.1 COUPLED TO FERMIONS AT FINITE DENSITY 

We can couple the Xcj)^ model to fermions yielding a mechanism for placing the system in a finite fermion 
density environment. This procedure may produce another mechanism for the symmetry restoration. In the 
dynamical picture we suppose that the coupling of fermions to the scalar field is switched off at t = as a first 
approximation to the problem. In part, this can be justified in the case one has decreasing scalars- fermions 
interaction amplitude with relation to the self interaction of the scalar field. If we consider an expanding 
finite density environment ("fireball") the density is expected to decrease as the system expands. As we will 
not be concerned with the fermion dynamics we will only consider the interacting part of its Hamiltonian. 

The following Hamiltonian density terms for spin half fermions '0 coupled to the scalar field are consid- 
ered: 

Hi = 5a0(x)V^(x)V(x) + (27) 

where the coupling constants ga and are dimensionless only in 3+1 and 2+1 dimensions respectively. 
Nevertheless we are allowed to consider them as effective couplings of an effective theory. Eventually one 
would need other couplings as one considers higher energy processes. 

The above couplings lead to changes in the equations of the Xcp^ model. In particular, we are interested 
in possible effects in the minimum energy state of the model. For this, we repeat the calculation of the GAP 
equation which is obtained by the minimization of the Hamiltonian density (expression |5|) in the frame of 
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the Gaussian approximation. Considering that the wave functional of the system now acquires a fermionic 
sector |H > (which may be a Slater determinant and is a function of the chemical potential) due to the 
presence of fermions at finite density we write: 



^[(p] >^ \^[(p] > x\E[i;]> . 



(28) 



With the fermionic wave functional one calculates the averaged values of the interacting Hamiltonian terms 
which enter in the effective potential of the scalar field and modify the GAP equations. These averaged 
terms are given in terms of a fermionic density: 



The density pj can be calculated by fixing the chemical potential eventually at a given temperature. However 
it is not the aim of this work to perform such a detailed dynamical self consistent microscopic description. 
Here we only want to provide a physical basis for the initial conditions of the time dependent model. By 
the minimization of the energy density with respect to the one and two point functions of the scalar field, 
we obtain the GAP equations (at finite fermionic density) which can be written as: 



Regularization of the ultra violet divergences and the subsequent renormalization of the parameters of the 
model does not change qualitatively these equations. 

We fix a coupling constant (a strong one: A ~ where fj, is the physical mass) for an one dimensional 

system but the conclusions do not change qualitatively for higher dimensions. Solutions for these two 
expressions (for fixed couplings) have been searched and the main conclusions due to the introduction of 
the fermionic density dependence through the couplings ga and gb are the following: 

(i) Keeping gb = and a quite small ga = +/i/5. As can be seen in expression ( |30| ) there will not exist 
the symmetric solution c^o = ^-t non zero density. In fact as the density increases the (real) condensate 
value will increase. This leads to a still more asymmetric phase. 

(ii) On the other hand for the same value of the other coupling gb = +^/5 but ga = the asymmetric 
phase will disappear when p ~ 14//^ and there will be no more condensates. At nearly double densities the 
physical mass also disappears. 




(29) 




(30) 
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(iii) If one assumes negative coupling gt < (and Qa = 0). In this case the condensate has increasing 
values for higher densities until a point where there is a complete disappearance of the asymmetric phase. 
This point coincides with a zero value for the mass from the GAP equation at p ~ 25 fj?. 

For finite temperature field theories, as 0(N) or 0(N)XO(N) models, the spontaneously symmetry 
breakdown may be restored or not at high temperatures ||3^. We point out, however, the possibility of 
further symmetry breaking at finite densities. A more complete analysis the finite density effects on scalar 
models will be shown elsewhere [^]. 

We obtained, therefore, another mechanism for considering local variations of the condensate: (p ^ (f)Q. 
This allows for scenarios in which one obtains non homogeneous configurations for the scalar fields at finite 
fermionic densities and temperatures fixing the initial conditions. The situation in which the condensate is 
suppressed is usually more accepted. We do not neglect the possibility of enhancing the condensate due to 
some different phenomena when the field interact with matter. 

4.2 CHOICE OF PARAMETERS 

We consider static initial conditions, the initial "velocities" of the classical and quantum parts of the field 
are taken to be zero in all examples shown below: 

SocG' = 7roc(^ = 0, (31) 

where the dot means time derivative. 

The model has two parameters which must be fixed: mass and coupling constant. We have chosen 
some values to place the system in the scaling limit. With the lattice spacing Ax = 0.1 fm we had considered 
Ax << ^ < L where ^ = 1/^ is the correlation length and L is the size of the lattice. In this region the 
universal properties of the lattice model can, in principle, be described by a continuum field theory. As a 
rule, the physical mass was chosen to be ;U = WOMeV for the dynamical situations. 

For the coupling constant, which has dimension fm""^ in one dimensional space, different values were 
chosen: from A = l/12/i^ ~ 0.021/m~^ to A = 600/m-~^ ~ 2350;U^. For low dimensions the model is 
super-renormalizable and the coupling constant is large [^ ]. The smaller value was already considered to 
be in the non perturbative regime [^]. We have found, however, that this is not the case for the examples 
shown in the present work. Indeed, we have found that the dynamics for couplings with values up to 
A ~ Ifm? ~ 5/12;U^ are not substantially different from the tree level case in the early times dynamics. 
Numerical examples will be shown in section 5. 
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In order to perform consistent comparisons between classical and quantum dynamics as well as among 
different initial conditions it is important to fix the parameters of the model. But instead of fixing physical 
mass and coupling, it is also possible to fix other variables as the energy density or particle number. For 
the sake of the argument, let us consider the simpler case of homogeneous solutions in the vacuum. At the 
tree level {G = 0) we obtain, in the vacuum: 

W.ac = -^. (32) 

Therefore we can fix, for example, the mass and the energy density of the vacuum and calculate the corre- 
sponding coupling constant. This can be useful for the study of the influence of the quantum effects because 
the inclusion of fluctuations (perturbatively or not) change the ground state which is deflned by the mass 
and coupling constant [jl^ ]. 

Let us consider a regularized energy density TCcvac at the Gaussian level which can be particularly well 
suited for the lattice calculations. Fixing the physical mass we obtain G. Writing the total energy density 
with expressions (^) and ( |Tl| ) we obtain a second degree polynomial expression for the coupling constant. In 
terms of the bare (regularized) quantities the corresponding solutions for the coupling constant as a function 
of the (regularized) energy density are given by: 



A = (33) 

where S = STLvac — G~^ + 2fi'^G. The renormalized version of H (as derived, for example, in |3^]) can also 
be used for this calculation. In one dimension there may have negative and complex solutions which do 



not seem to correspond to meaningful stable minima of the effective potential |16|. This calculation could 
be meaningful even in 3+1 dimensions to the extent that the Xcp^ model can be considered as effective for 
which the cutoff can be fixed at some high energy scale. In this case the bare and renormalized quantities 
can be related by expressions (|l^ ) in one spatial dimension. This is equivalent to placing the system in a 
lattice, which provides us with a natural cutoff (l/Ax), and to perform all calculations on it. 

4.3 THE PICTURE 

We have therefore the following picture. Firstly we fix two parameters of the model, which will allow us 
to make meaningful comparisons between the (time dependent) tree and Gaussian levels. Secondly, it is 
assumed that the scalar field locally experiences an interaction with a finite fermion density in a small region 
or has some contact with a thermal bath. These interactions -which change the ground state of the model- 
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are switched off at t=0 yielding non iiomogeneous (and non equiiibrium) initial conditions for tlie scalar 
field which are evolved in time. The temporal evolution is performed within the tree level and Gaussian 
approach equations, producing the expansion of regions (bubbles) endowed with high energy densities. 



4.4 INITIAL CONDITIONS: FINITE DENSITY AND TEMPERATURE 

Firstly we suppose that the temperature varies continuously from the center of the lattice, where there is 
a high energy region, to the (zero temperature) vacuum. This constitutes an out of equilibrium situation 
which can be implemented in a lattice by the following configuration: 

^'^(x^t = 0)= fj.hanh' (^^^) • (34) 

The bubble of high energy density is centered at x = L/2 (L being the size of the lattice) and it has size 
given by A (which will be taken to be 0.5/m). At the center of this bubble = OMeV which implies 
very high temperatures. As the temperature also modifies the order parameter of the model {(p) it may 
eventually yield a symmetry restoration. 

We will consider another kind of configuration given by: 

<^(x, t = 0) = ^otanh^ f ^— • (35) 



A 

In the central region of the lattice the condensate is suppressed. This can happen due to the interaction 
with a finite density medium. An unifying picture can be associated to the above configurations. In high 
energy collisions there is a large amount of energy deposited in a small region, eventually making the system 
to have its parameters modified, such as particle masses and condensate. This energy excess is expected to 
propagate (expand) in the space time. 

An enhancement of the condensate due to the interaction with matter can also be assumed (it could 
correspond to further symmetry breaking at high densities): 

2 



(x,t = 0) = (po 



1 + asech 



(36) 



A 

Where a is a real positive number which measures the amount of energy excess deposited in the central 
region of the lattice. 

Finally, another kind of initial conditions will be considered. It corresponds to a deviation from the 
KINK classical solution to the equations of movement. The static KINK solution of the classical A</)^ model 



is given by |35]: 



(x, t = 0) = ^otanh ( ) , (37) 
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where B = y—2/mQ. In this case cj){r = 0) = — 0o ^-nd cj){r = 10) = c^O) where r = and r = 10 are the 
borders of the lattice. This configuration can be seen as a wall which separates regions with different vacua 
where (j) = i^o- It is stable at the tree level and corresponds to a non trivial topology case. By considering 
anti-periodic boundary conditions we have chosen a deviation with relation to the kink solution: 

<^(x, t = 0) = 4)otanh (^—^^^ ■ (38) 

Where B was considered to be i? ~ 1/(4^). There is an energy excess in the central region of the lattice, 
over the Kink. 



5 NUMERICAL RESULTS 

In Figure 1 some of the initial conditions for the field {G{r) and 4){r)) are shown. In thick solid line it is 
shown G = G{m?{r),t = 0) corresponding to a bubble of high energy density which comes out from a zero 
mass region inside the vacuum (for // = ISOMe^). The physical mass is given by expression 



fx — L l2\ 

rn^it = 0) = fi^tanh^ i — 05 j ' (t>{t = 0) = (po- 



In this case, it is therefore assumed that the deviation of the mass in the center of the lattice but the 
condensate is kept constant. The other two curves (thin solid and dotted lines) correspond to non trivial 
initial conditions for the condensate ( |35| , p6[ ) and will be discussed later. 

The first initial condition ( p^ is evolved for two different values of the coupling constant: A = /i^/12 ~ 
0.02/m~^ and 600/m"^, for a physical mass /x = lOOMe^. The resulting evolution of the energy density 
configuration for low times are shown respectively in Figures 2 and 3. In figure 2 the energy density 
configuration is exhibited at different time steps: t = 0.1/m, 1.0/m, 2.0/m and 3.5/m as a function of 
the lattice points. The initial (higher) energy bubble expands (by "waves") towards the extremities of the 
lattice distributing the initial "potential energy" over the lattice and among classical and quantum degrees 
of freedom. For a much stronger coupling constant, figure 3, the energy density configurations are shown 
and expand in a more "concentrated way" than in the weaker coupling case. Almost no difference can be 
seen in what concerns the expansion velocities. 

The following initial condition is considered for the next cases (expression ([35|)): 

0(t = 0) = ^otanh^ (^J^) ' ^ = ^- 
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In this case the condensate is suppressed in the central region of the discretized space-time. In figure 4 the 
evolution of the corresponding energy density is shown for a quite weak coupling constant A = 0.6/m~^ ~ 
2.36/i^ without quantum fluctuations, i.e., by means of the classical equation of motion -third and fourth 
equations of (0). Two-peak waves expand from the center of the lattice to the borders. The initial energy 
density has a two-peak structure which happens due to the fact that the gradient term in the Hamiltonian 
density has the largest values (see expression (^)). This figure can be compared to the evolution of the 
complete system including quantum fluctuations which is shown in figure 5. The two point function G is 
considered to be in the vacuum value at t = 0. Due to the weakness of the coupling constant (the same as 
the one used in figure 4) there is no visible difference between both figures. The main difference, which is 
quite big but not relevant to the dynamics, is the overall normalization of the energy density. The inclusion 
of quantum fluctuations introduces a zero point energy whose effect on the dynamics is not relevant, as 
the comparison between figures 4 and 5 shows. The above value of the coupling constant still is in the 
(dynamical) perturbative regime for this early time analysis. Comparison of this weak coupling constant 
case (A ~ 1/12;U^) to other works [11 1 suggests that the non perturbative dynamics may be present for some 



kind of configurations but not for others in the early time analysis. 

We consider now another type of initial condition. The value of the condensate is enhanced in the central 
region of the lattice by means of expression ( |36[ ) : 

a fx - 



ct){t = 0) = </.o 



1 + asech 



A 



where a is a real number which will be assumed to be in the range < a < 0.5. It measures the energy 
of the bubble which is placed in the center of the lattice. It still has size roughly given by j4 = 0.5/m. 
Since the previous examples have showed that A = ~ is not strong enough as to yield non 

perturbative quantum effects we consider a stronger value (A = ~ 236^^) for the next figures. This 

configuration could correspond either to a case the model is place in a region with finite fermion density 
with the particular kind of coupling discussed in the previous section or to a situation where external fields 
would "merge" in a condensate in that small region increasing its value and the respective energy density 

Hi- 

In figure 6 the classical energy density is shown for an initial condition given above with a = 0.2, which 
represents a not very energetic bubble. The energy density evolves nearly as the same way as the case of 
figure 4 (initial condition defined by a suppression of the condensate instead of enhancement). There is 
a small difference which is related to the relative amount of energy in each peak of the expanding waves 
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at larger times (when t ~ 3.5/m). In the present case (enhancement of the condensate) the second peak 
is shghtly higher with relation to the configuration where the condensate is suppressed (figure 4). By 
switching on quantum fluctuations we obtain the temporal evolution shown in flgure 7. The first issue we 
can note is that between time steps 0.1/m and 1.0/m there is an "amplification" of the locally concentrated 
energy density leading regions in which its value is smaller than the vacuum energy density. Theses regions 
expand. These are local effects since the total energy is always conserved and given by the vacuum energy 
plus the initial bubble energy excess (positive). Furthermore, these regions tend to disappear for larger 
times. Another important effect is that the most part of the expanding energy is now concentrated in the 
"advanced" peaks, which arrive first in the extremes of the lattice. This corresponds to an acceleration of 
the expansion |jl5|. For the preceding cases the energy density "bumps" took nearly At ~ 5 fm to arrive at 
the borders of the lattice. In the present case (figure 7), this time is reduced to nearly At ~ 3.8 fm. As noted 
in this last reference the energy density becomes smoother when considering the quantum fiuctuations but 
this happens only in the non perturbative region of the coupling constant. It is also possible to note that 
the central region of the lattice, from where the bubble started expanding, tends to assume energy density 
values close to that of the vacuum as the energy goes away towards the borders of the lattice. Stronger 
coupling constants therefore lead to the following quantum effect: there is an acceleration of the energy 
density expansion. Although the two-peaked form of the expanding waves persists, the advanced peaks 
(those which arrive before to the borders of the lattice) are strengthened with relation to the others. This 
is a clear indication that the expansion is faster corresponding to a non perturbative quantum effect. 

What happen if the initial energy density excess in the center of the lattice is increased? In figures 8 
and 9 this is shown respectively for the classical system and for the "complete" (classical plus quantum) 
field with the same set of parameters as the preceding figures but considering a = 0.5 in expression (^6|). 
The coupling constant is in the range A = 60/m-2 ~ 236/i2. The energy density excess is much higher 
than those considered before as it can be checked by comparing with figure 8. In this figure the temporal 
evolution of the energy density of the condensate without quantum fiuctuations is shown. In spite of different 
normalization and energy values the classical dynamics is not modified with relation to figure 4 (for which 
the initial condition was given by the suppression of the condensate with expression ( |35[ ) ). In figure 9 the 
dynamics of quantum fluctuations is taken into account. It is also noted in this case a quite large energy 
density amplification in the beginning of the evolution and the existence of expanding regions in which the 
energy density is smaller than the value of the vacuum. They modify the two-peak structure discussed 
before. Although the last step in such time evolution present values for the energy density which may not 
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be completely reliable due to numerical uncertainties the main issues are completely trustful even because, 
we emphasize, the total energy is conserved. The acceleration found before still is present. Moreover, it 
is interesting to note that there is little difference between the temporal evolution of cases in which initial 
conditions are given in term of a with different values, i.e., for a = 0.2 or a = 0.5 in the early time dynamics. 

For the examples shown above only the short time behavior was analyzed. These are so short time scales 
that the boundary conditions are not even relevant. The large time behavior was briefly studied for some 
cases in the strong coupling limit (A = 60/m~'^). No equilibration was found for larger times evolutions, 
i.e., the amplitude of the (classical and quantum) oscillations do not tend to zero. There are several other 
different approaches dealing with different aspects of the non-equilibrium field dynamics 21, 3f:]. 

In the next figures another kind of (non trivial) initial configuration is considered for the condensate. 
In this case <^(r = 0) = — 0o and 0(r = 10) = where and 10 are the borders of the lattice. This 
configuration can be seen as a wall which separates regions with different vacua where = ±i^o- By 
considering anti-periodic boundary conditions we have chosen a deviation with relation to the KINK solution 
(expression |^) p5| ]: 

^{x^ i = 0) = ^Qtanh • (39) 

Where B was considered to be i3 ~ l/(4/i) and A = 60/m~^. Firstly the evolution of the equations 
of movement for the condensate without quantum fiuctuations (G = 0) is performed. In figure 10a the 
resulting classical field profile {(f>{t)) at the points x = bfm and x = 10/m are shown. The change of the 
field at (j){x = 10, t) (there is a "flip") is due to the anti-periodic boundary conditions and happens exactly 
at the time when the energy density bumps arrive at the border of the lattice. The field in the central 
region has a static and constant value given at t = 0. In figure 10b we show the energy density expansion (it 
does not have the two-peak structure due to the initial condition). By switching on quantum fiuctuations 
we obtain figures 11a, b, c. In figure 11a the condensate in the same points of the lattice as in figure 10a 
is shown. The dynamics is similar but the field (at r=10 fm) "flips" to the other vacuum value faster. In 
figure lib the deviation of the quantum fiuctuations with relation to the value in the vacuum, at the same 
points {6G{t,r) = G{t,r) — Go), are shown. They clearly exhibit the energy transfer between classical and 
quantum degrees of freedom. In particular, when the "condensate fiip" occurs, the quantum fiuctuations 
are enhanced in the corresponding point (r = 10 fm). For the energy density configurations, in figure 11c, 
there is no strong effect with relation to the purely classical dynamics (seen in figure 10b). 

For the figures 12 and 13 a different size of the bubble of suppressed condensate is considered with 
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the initial condition given by expression (^). We have considered half value, i.e., A = 0.25/m. This 
yields a smaller region out of equilibrium. In figure 12 the temporal evolution of the energy density of the 
classical condensate case is shown as well as in figure 13 the same for the classical plus quantum system. 
The concentration of energy density in the "advanced" peaks discussed before (figures 5,7 and ||l^) is well 
visible in figure 13. This indicates, again, an acceleralation of the expansion with relation to the classical 
dynamics. Moreover, there is the appearance of new bumps in the expansion due to the inclusion of the 
quantum dynamics in a finite (discrete) system. 

In the continuum limit the GAP equations and the equations of movement present the ultraviolet (UV) 
divergences discussed in section 2. This is reproduced in the lattice calculations and to show that our 
results do reproduce features of the continuum renormalized model we have studied the limit of smaller 
lattice spacings down to Ax = 0.02/m. In order to continue in the scaling limit the mass /i^ must be 
changed in the same way as the GAP equations in the lattice since it absorbs the UV infinities in the 



renormalization. The two point function scales as |33]: 

where is the field renormalization factor (it is finite in the 1+1 dimensional case). Results for the small 
lattice spacing limit are not visibly modified. The difference is found in the normalizations (absolute values) 
of the two point function and of the energy density. By subtracting these values by the vacuum ones the 
results remain unchanged. Besides that, the dynamics is not affected by these divergences, i.e., by the 
regularization method. This was shown in section 2.3 for the particular case of initial conditions given by 
small deviations from the vacuum. 



6 SUMMARY 

We have analyzed the temporal evolution of expanding non-homogeneous configurations of the Xcj)'^ model 
considering two different approaches: the classical equations of motion and compared its results to the 
equations of motion in the frame of the Gaussian variational approximation in a 1+1 dimensional lattice. A 
schematic model for the model at finite fermionic density has been drawn for the initial conditions and the 
equations of movement have been solved. The condensate may either disappear (symmetry restoration) or 
become higher (no symmetric phase, with further symmetry breaking) at higher densities. The parameters of 
the A(/)^ model were fixed in order to allow a comparison between the classical and quantum field temporal 
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evolution. We have been able to study the influence of the quantum fluctuations on the classical field 
dynamics for different sets of free parameters. By varying the parameters of the model and the non- 
homogeneous initial conditions we have carefully investigated the expansion of different field configurations 
in the frame of the Gaussian approximation. The quantum fluctuations accelerate the expansion of a 
concentrated configuration of the field. This effect is considerable for strong coupling constants and particular 
cases of the initial conditions, namely when there is an enhancement of the condensate (stronger symmetry 
breaking) instead of suppression (restoration). However, no big differences were found for these two different 
initial conditions in the early time dynamics. Closely related works have been done by 11, respectively 



for other initial conditions, as for instance a Gaussian configuration, and additional averaging over ensembles 
of mean fields (which seems to lead to thermalization at long times). It is possible to conclude that the 
initial conditions play an important role in the temporal evolution. We have been concerned mainly with 
short time intervals evolution and other issues related to thermalization and equilibration have not been 
addressed extensively. 
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Figure Captions 

Fig. 1 Examples of initial field configuration. Thick solid line corresponds to the two point function G(r, t = 0) 
for the initial condition ( |3^ considering fi = 150MeV. Thin solid and dotted lines correspond to the 
initial condensate configuration of expressions ( |35[ ) and (36 with a = 0.2) respectively. 



Fig. 2 Evolution of the energy density of the field (classical + quantum) at different times. The initial time 
mass configuration is given by (^) and A = ^^/12. 

Fig. 3 The same as Fig.2 with A = 600/m~^. 

Fig. 4 Evolution of the energy distribution of the classical field for the initial configuration (^) and A = 
0.6/m~2. 
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Fig. 5 Evolution of the energy distribution of the field with the quantum fluctuations for initial configuration 
dD and A = 0.6fm~'^. 

Fig. 6 Evolution of the energy distribution of the classical field for the initial configuration (^) with a = 0.2 
and for A = 60/m~^. 

Fig. 7 Evolution of the energy distribution of the field with the quantum fluctuations for initial configuration 
(H) with a = 0.2 and A = 60/m-2. 

Fig. 8 Evolution of the energy distribution of the classical field for the initial configuration (|3^ ) a = 0.5. 

Fig. 9 Evolution of the energy distribution of the field with the quantum fluctuations for initial configuration 
dD with a = 0.5 and A = 60/m~^ 

Fig. 10 Evolution of the energy distribution of the classical field (f){t, x) for the initial configuration (|37| ) and 
A = 60/m-2. 

Fig. 11 Evolution of the energy distribution of the field with the quantum fluctuations for initial configuration 
@) with A = 60/m-2. 

Fig. 12 Evolution of the energy distribution of the classical field for the initial configuration (|35| ) with A = 
0.25/m. 

Fig. 13 Evolution of the energy distribution of the field with the quantum fluctuations for initial configuration 
(H) with A = 0.25/m. 
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